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(57) Abstract 

The present invention relates to a method 
and system for mobile location estimation in which 
range measurements between the mobile station 
and a base station are discriminated as being ei- 
ther from a base station which is line of sight 
with the mobile station or a base station that is 
non-line of sight with the mobile station (22). 
Non-line of sight error present in range measure- 
ments from base stations that are non-line of sight 
with the mobile station is corrected (24). A recon- 
structed line of sight range measurement is estab- 
lished from the correction of the non-line of sight 
error. The reconstructed line of sight range mea- 
surements can be used with range measurements 
of base stations determined to be line of sight for 
accurately locating the mobile station (26). 
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METHOD AND SYSTEM FOR MOBILE LOCATION ESTIMATION 



This application claims the benefit of U S Provisional Application No. 60/027,453 
entitled Non-Line Of Sight Problem in Mobile Location Estimation filed by Applicants on 
September 27, 1996 hereby incorporated by reference into this application 

Background of the Invention 

1 Field of the Invention 

The present invention relates to a method and system for mobile station location 
estimation in which base stations that are in line of sight of the mobile station and base stations 
that are not in the line of sight of the base station can be determined. Errors in base station 
signals generated from determined non-line of sight base stations are reduced for providing 
improved mobile station location estimation. 
2. Description of the Related Art 

Mobile location estimation determines a geographical estimate of the location of a 
mobile station Mobile location estimation is useful in management of fleets of mobile 
stations, location dependent information services, location dependent billing services and 
Emergency 911 location of a mobile station Enhanced 911 is designed to automatically 
forward the number of a caller to a public safety answering point (PSAP) In implementing 
enhanced 911 in a wireless network, wireless service providers provide two dimensional 
location of the vehicle to the public safety answering point (PSAP) The Federal 
Communications Commission (FCC) has regulated by the year 2001 that wireless service 
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5 providers have the capability of locat.ng callers in two dimensions within 125 meters 67% of 
the time. 

One conventional method for locating a mobile station in two dimensions would use 
the measurement of the line of sight distance between the mobile station and at least three base 
stations U.S. Patent No. 5,365.516 describes a method for determining the location of a 

10 transponder unit in which a radio signal is sent by the mobile station The arrival time of the 
rad.o s.gnal is measured at each of three base stations. Each distance measurement between 
the mobile station and one of the base stations can be used to generate a circle wh.ch is 
centered at the measuring base station. The circle has a radius which is equal to the distance 
between the mobile station and the base station Accordingly, three circles are generated, one 

1 5 for each of the base stations. In the absence of any measurement error of the distance between 
the base stations and the mobile station, the intersect.on of the three c.rcles unambiguously 
determines the location of the mobile station This method has the drawback that the distance 
measurements can be corrupted by noise resulting in errors in determining the locat.on of the 
mobile station. 

20 A conventional solution for providing more accurate position estimates is to reduce the 

error due to noise with a least squares analysis Accordingly, the least squares analysis 
provides a more accurate position estimate This solution has the limitation of not accounting 
for the possibility of a lack of a direct path between the base station and the mobile station 
For example, in an urban environment, a building or buildings may be in the path between the 

25 mobile station and the base station A propagating signal between the mobile station and the 
base station can be reflected and defracted by the object in the path of the mobile station to the 
base station resulting in the signal traveling excess path lengths The excess path lengths can 
be on the order of a hundred meters. 
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5 The lack of direct path between the mobile station and the base station can be defined 

as a non-line of sight (NLOS) The importance of detecting and reducing the NLOS 
measurements between a mobile station and a base station is recognized in M.I Silventoinen. 
et al., "Mobile Station Locating in GSM', IEEE Wireless Cojmimnication System Symp osium 
Long Island NY, November 1995 and J.L. Caffrey et al., "Radio Location in Urban CDMA 
10 Microcells", Proceedings of the Personal In door and Mohile Radio Enyjnmment 1 QQ^ 

U.S. Patent No 5,365.516 ('516 patent) describes an embodiment of a transreceiver 
locating system operating in an environment susceptible to multipath interference. The system 
includes a transponder which is operable within a prescribed coverage area to transmit a burst 
of data symbols in a coded carrier pulse. Each base station includes a receiver for detecting 
15 and responding to the data symbol at a given time, interrupting the data symbol and rejecting 
echoes resulting from multipath interference. A comparison circuit responds to the receiver 
for comparing respectively identified given times and decorrelating the t.me difference to 
improve data quality Although the "516 patent addresses multipath .nterference. it does not 
attempt to detect base stations for reducing multipath NLOS with mobile stations 
20 It is desirable to provide a method and system for providing improved mobile location 

estimation which is robust to NLOS error. 

Summary of the Invention 
Briefly described, the present invention relates to a method and a system for mobile 
location estimation in which base stations are identified to be either line of sight (LOS) or non- 
25 line of sight (NLOS) with a mobile station. A range measurement is determined as the distance 
between the base station and the mobile station. NLOS ranging error is corrected for base 
stations identified to be NLOS with the mobile station by reconstructing the LOS 
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5 measurement. From the range measurements of base stat.ons ident.fied as LOS and the 
reconstructed LOS range measurements the location of the mobile station is estimated. 

The base station can be identified as being NLOS by comparing the standard deviation 
of standard measurement noise from the environment to the standard deviation of a smoothed 
range measurement obtained from the range measurements between the base station and 

10 mobile station. The smoothed range measurement can be obtained using an N* order 
polynomial fit. It has been found that when the standard deviation of the smoothed range 
measurement is on the order of the standard deviation of the standard measurement no.se, the 
base station corresponds to an LOS env.ronment and when the standard deviation of the 
smoothed range measurement is greater than the standard deviation due to standard 

1 5 measurement noise, the base station corresponds to an NLOS environment Alternatively, the 
residuals from a least squares analysis can be used to determine the presence of NLOS range 
measurements. 

NLOS error can be corrected when the standard measurement noise dominates the 
NLOS error and there is predetermined identification of the approximate support of the 

20 standard measurement noise over the real axis. A reconstructed LOS range measurement can 
be determined by graphing a curve of the smoothed range measurements. The point of 
maximum deviation of the smoothed range measurement below the curve is determined. The 
curve is displaced downwards to pass through the point of maximum deviation Thereafter, 
the curve is displaced upwards by the value of the maximum standard measurement noise 

25 deviation from an LOS measurement with negligible noise, thereby providing a reconstructed 
range measurement. 

The mobile location estimation can be determined using at least three range 
measurements between LOS base stations and the mobile station or reconstructed LOS range 
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5 measurements in a muitilateration analysis. In this analysis, a circle is generated from each 
range measurement The circle is centered at the base station and the range measurement is 
the radius of the circle. The estimated intersection of the three circles determines the location 
of the mobile station. Alternatively, two range measurements and information directed to the 
position angle of the mobile station can be used for estimating the location of the mobile 
1 0 station. 

The present invention has the advantages of accurately determining the location of a 
mobile station by reducing NLOS error. In addition, the present invention can provide 
confidence in an LOS environment that all base stations are LOS with the mobile station 
Results indicate that position range bias due to NLOS error can be reduced several orders of 
15 magnitude with the method of the present invention. 

The present invention will be more fully described by reference to the following 
drawings. 

Brief Description of the Drawings 
20 F| S 1A is a schematic diagram of an environment in which there is an unobstructed 

line of sight radio signal path between a mobile station and a base station 

Fig. IB is a schematic diagram of an environment in which there is a non-line sight 
radio signal path between a mobile station and a base station 

Fig. 2 is a flow diagram of the system and method for mobile location estimation in 
25 accordance with the teachings of the present invention 

Fig. 3 is a schematic diagram of distance measurements of a reconstructed line of sight 
base station and determined line of sight base stations. 
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5 Fig. 4 is a flow diagram of a method for identifying non-line of sight base stations of 

the present invention. 

Fig. 5 is a flow diagram of an alternate method for identifying non-line of sight base 
stations. 

Fig. 6 is a flow diagram of a method for reconstructing a line of sight base station for 
1 0 non-line of sight measurements. 

Fig. 7 is a graph of a comparison of NLOS measurements and reconstructed LOS 
measurements. 

Fig. 8 is a schematic diagram of a system for implementing the method of the present 
invention 

15 Fig 9 is a schematic diagram of positioning of base stations used in examples of 

performance of the method of the present invention 

Fig. 1 OA is a graph of two dimensional tracking without non line of sight error 
detection and correction 

Fig. 10B is a graph of two dimensional tracking with non line of sight error detect.on 
20 and correction. 

Fig. 11A is a graph of two dimensional tracking without non line of sigh: error 
detection and correction 

Fig. 1 IB is a graph of two dimensional tracking with non line of sight error detect.on 
and correction. 

25 Fig. 12 is a graph of the fraction of time a base station was declared NLOS ,sing the 

residual rank analysis method 
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5 Detailed Description of the Present Invention 

During the course of this description like numbers will be used to identify like elements 
according to the different figures which illustrate the invention. 

Fig. 1A illustrates a schematic diagram of a line of sight (LOS) path 10 between a base 
station 12 and mobile station 14. Signal 13 can be transmitted from base station 12 to mobile 
10 station 14 and returned from mobile station 14 to base station 12 Fig IB illustrates a 
schematic non-line of sight (NLOS) path 1 1 between base station 12 and mobile stat.on 14 
Building 15 is positioned between base station 12 and mobile station 14 resulting in reflection 
of signal 16 For example, signal 13 and signal 16 can be a radio signal. 

A range measurement for measuring the distance between base station 12 and mobile 
15 station 14 can be measured as the time it takes a signal sent between base station 12 and 
mobile station 14: 

r = cT (|) 
in which the mobile station to base station range measurement is represented by r, c represents 
the speed of light which is the same speed as the propagation of radio waves and T represents 
20 the one-way travel time of the signal A range measurement of the distance between mobile 
station 14 and base station 12 in Fig. 1 A and IB can be determined using equation (1) based 
on travel time of signal 13 and signal 16, respectively, between base station 12 and mobile 
station 14 The value of r generated from signal 16 is greater than the value of r generated for 
signal 13. 

25 Fig. 2 is a flow diagram of the system and metHod of the present invention for mobile 

location estimation 20. In block 21, a range measurement is obtained between mobile stat.on 
14 and base station 12 using equation (1) In block 22, base stat.on 12 is identified as being in 
line of sight (LOS) or non-line of sight (NLOS) with mobile station 14 Block 22 is repeated 
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for a plurality of base stations 12 positioned at different locations from mobile station 14 If 
base station 12 is identified to be LOS in block 22, the range measurement obtained from 
block 21 is forwarded to block 26 If base station 12 is identified to be NLOS in block 22. 
block 24 is implemented for reducing the error of the range measurement between base station 
12 and mobile station 14, thereby rendering the range measurement between base station 12 
and mobile station 14 as a reconstructed LOS base station 13, as shown in Fig 3. 

In Fig. 3, base station 12 labeled BS1 has a range measurement labeled LOS RANGE 
1 determined to be LOS Base station 12 labeled BS2 has a range measurement labeled LOS 
RANGE 2 determined to be LOS. Base station 13 labeled BS3 has a range measurement 
labeled NLOS RANGE 3 determined to be NLOS A range measurement for the 
reconstructed LOS base station labeled RECONSTRUCTED RANGE 3 is forwarded to block 
26 Range measurements from determined LOS base stations from block 22 labeled LOS 
RANGE 1, and LOS RANGE 2 are also forwarded to block 26 From the range 
measurements of the reconstructed LOS base stations or the determined LOS base stations, or 
a combination of range measurements of the reconstructed LOS the mobile location estimation 
can be identified using a conventional multilateration technique, such as described in U.S. 
Patent No. 5,365,516, hereby incorporated by reference into this application Alternatively, 
the mobile location estimation can be determined from time difference of arrival time 
measurements as the difference of propagation delays between the mobile station 14 and pairs 
of base stations 12. In this case, the position estimate is at the intersection of hyperbolas. The 
number of base stations can be reduced below three if there is also angle of arrival information 
These methods are described in IS. Rappaport et a!., " Position Location Using Wireless 
Communication On Highways Of the Future", IEEE Communication. M a f , a7 in„ October 
1996. 
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One method for identifying if a base station is LOS or NLOS in block 22 is illustrated 
in Fig 4. In this method, the time history of range measurements between base station 12 and 
mobile station 14 is combined with predetermined standard deviation from conventional 
measurement noise in a radio signal environment. 

The arrival time of signals sent from base station 12 to a mobile station 14 and 
transponded back to the base station 12 can be converted to a range measurement, in block 
30. The range measurement at the m" 1 base station at time t k can be represented as 
r„(t.) =L m {t.) + //„,(/,) + NLOS m (l.) (2) 

for /77 = 1, ...,M i = 0, ... K-l, wherein 

L m (/, ) is the LOS distance between a mobile station and the m* base station in two 
dimensions which is given by: 

L m {t,)= \x{t,)+j *y(t,) - Xm -j *y m \, (3) 
*('. ), >>(t, ) and (x m , y m ) are respectively the coordinates of the mobile station at time, t, , and 
those of the m' h base station; n m (t, ) represents conventional measurement noise such as 
additive white Gaussian measurement noise and NLOS„ (/,) represents NLOS measurement 
error at time t, .; and M is the total number of base stations; and A' is the total number of time 
samples. 

In block 30, an LOS range measurement with negligible noise is obtained for base 
station 12 in LOS with mobile station 14 The LOS range measurement can be obtained by 
physically measuring a range between base station 12 and mobile station 14 or can be obtained 
as a range measurement determined by equation (I) in a negligible noise environment. In 
block 31, a noisy range measurement is determined as a range measurement which is LOS 
with a base station taken in a noisy environment In block 32, the standard deviation of the 
noisy range measurement from the LOS measurement without noise is determined. Blocks 30, 
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31 and 32 can be predetermined before identifying base station 12 as either LOS or NLOS in 
block 22 The standard deviation due to noise n m (t) can be represented by a m . 

In block 32, the range measurement obtained from block 21 is smoothed by modeling 

A/-/ 

, r m (l.) = 2X(»K" (4) 



n~0 



and solving for the unknown coefficients, {a m (n)}":' with a least squares technique The 
smoothed range measurement can be represented as: 



n-0 



(5) 



In block 34, the standard deviation of the smoothed range measurement from a noisy 
range measurement (i.e., the residual) is determined The standard deviation of the residual 
from block 34 can be represented as 6 m since a' m = Efn' m </)). The smoothed range 
measurements along with the noisy range measurement can be used to determine standard 
deviation a m with the formulation of: 



^ = j7Z(*.(0-%(0)' 



K£,\»mV,l-r m (i,)) (6) 

^ I =0 



From the value of the standard deviation, a m and the standard deviation o„, the range 
measurement can be determined as either the result of base station 12 being LOS or NLOS, in 
block 36 When the range measurement has NLOS error, the value of the standard dev.at.on 
a m is significantly larger than the value of the standard deviation a m Accordingly, range 
measurement for base station 12 that is NLOS with mobile station 14 is determined when the 
a m is greater than the standard deviation o m A range measurement of base station 12 that is 
LOS with mobile station 14 is determined when the standard deviation d m is on the order of 
the standard deviation o B1 
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Alternatively, a residual analysis ranking method can be used to identify a range 
measurement as being from a base station 13 NLOS with mobile station 14. Range 
measurements between mobile station 12 and base station 14 which have been obtained in 
block 21 are inputted to block 41 At each instance of time t,. estimated coordinates 
*«(',). of mobile station 14 are determined as least squares estimates in block 41 
The estimated coordinates 

x is(*, ).}'ls( 1 , )) are selected to minimize the formulation: 

■i/ 

F , ^(UO -MO)-' (7) 

m~ 1 

where L m (/, ) = | f (/, ) - x m + j. y(t, )-J*y m \. 

In block 41, a calculated range measurement is determined from the estimated 
coordinates. In block 42, a residual difference of the range measurement between mob.le 
station 12 and base station 14 with the calculated range measurement is determined The 
residual difference can be represented as: 

emit, ) = /"„,(/,)- L m (/, ) ( 8 ) 

In block 44, the number of times the residual difference of a range measurement to a 
base station 12 has the largest value in comparison to the residual difference determined for 
range measurements at other base stations is counted for each time instant /,. It has been 
found that base stations having a range measurement between a base station NLOS with a 
mobile station have a significantly larger number of greatest absolute residual differences than 
the number of greatest absolute residual differences from other base stations From the value 
of the counted number of residual differences, base station 14 can be defined as a base station 
12 that is a LOS or a base station 12 that is NLOS with mobile station 14 
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Fig. 6 illustrates a method for correction of range measurements between a base 
station 12 that has been determined NLOS with mobile station 14 to reconstruct a LOS range 
measurement. Data related to the range measurements from block 21 are smoothed using an 
N* order polynomial fit described in block 32. The smooth range measurements are inputted 
to block 52 The maximum deviation below the smoothed curve due to NLOS error in 
determined in block 56. It has been found that NLOS error is a non-negative random variable 
which can be approximately represented in a real axis as follows: 

0 <_ NLOSm (/, ) < p m 

in which Pm is the maximum value of NLOS error The standard measurement no.se, n m (t,) 
can be represented as a zero-mean random variable wh.ch can be approximately represented in 
a real axis as follows: - a m <;n m (t.) s a m , so that in a range measurement in which there is also 
an NLOS error, the total noise component can be approximated represented over the real axis 
as follows: 

-a. s n m (0 + NLOS m (/) <;0 m - Qm 

It has been found that the point of maximum deviation of the measured range below 
the smoothed curve is about a m below the LOS function represented as L m (t,). In block 58. 
the smoothed curve is displaced mathematically downward to the point of maximum deviat.on 
The smoothed curve is displaced mathematically upward by a value of the noise deviation a m 
in block 60 to provide a reconstructed curve representing a reconstructed LOS base station 

Fig. 7 represents a graph of a comparison of simulated range measurements Curve 90 
represents the true time range measurement between a' base station 12 which is LOS with a 
mobile station 14 Curve 91 represents determined range measurements having NLOS error 
Curve 92 represents a smoothed range measurement of block station 12 and mobile station 14 
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determined from block 30 of Fig. 4 Curve 93 represents base station 12 which is 
reconstructed LOS with the mobile station 14 from block 60 of Fig 6 

Fig 8 is a schematic diagram of a system 80 for implementing the method for mobile 
location estimation System 80 includes base station server 81 Base station server 81 can be a 
computer located at base station 12 or networked thereto Base station server 81 
communicates with base station 12 for requesting and receiving data related to range 
measurements of mobile station 14 and base station 12 Base station server 81 also collects 
information on range measurements between mobile station 14 and each of base stations 81 A- 
8 IN The information is reported to base station server 81 either by mobile station 14 or base 
station servers 81 A-8 IN. The functions of modules shown in Figs 4-6 which are coded with a 
standard programming language, such as C" programm.ng language. The coded modules can 
be executed by base station server 8 1 

Results for examples of mobile location estimates with system 80 are shown in Tables 1 
- IV and Fig. 9 through Fig. 12. In all of the examples, the veh.cle s position in the x-y plane 
at any is given by: 

x(t) = x„ + vj 

AO = }'o + v y i 
x(t) represents the x - coordinate in x -y plane at time instant, t, 
X') represents the y - coordinate in r - v plane at time instant, t, 
x„ represents the initial x - coordinate, 
y 0 represents the initial y - coordinate, 
v x represents the speed in x - direction, 
v y represents the speed my - direction 
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The sampling period was chosen to be 0.5s and 200 samples were taken. The velocity 
remained constant at v, - 9 7 m/s and v y = .6.8 m/s. Base stations 12 were assigned to have 
NLOS or LOS range measurements. Standard deviation of the standard measurement noise 
was represented as a a m was 1 50m and B m was chosen as 1 300m In each example three base 
stations 1 0 1 , 1 02, 1 03 were used uniformly spaced around a circle of 5 kilometers and a fourth 
base station 104 was located at the center of the circle, as shown in Fig. 9 

In a first example, base station 101 and base station 102 provide NLOS range 
measurements and base stations 103 and base station 104 prov.de LOS range measurements 
The standard deviation a m (m, of the smoothed curve determmed in Fig 4 is shown in Table 



TABLE 1 

STANDARD DEVIATION OF MEASUREMENTS FROM SMOOTHED 
CURVE FOR 2 NLOS MEASUREMENTS 



Base 


NLOS 


° m (m) 


101 


Yes 


467.3 


102 


Yes 


447.6 


103 


No 


163 1 


104 


No 


142.1 



The results indicate base stations 101 and 102 have NLOS range measurements with a 
significantly larger standard deviation than base station 103 and base 
range measurement. 



station 104 having a LOS 



Fig. 10A shows two-dimensional tracking error without NLOS identification and 
correct.on Fig. I0B shows two dimensional tracking error after the method of mobile location 
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estimation of the present invention is performed. The results indicate improvement of 
estimated vehicle trajectory after NLOS identification and correction. 

In a second example, base stations 101, 102, 103 and 104 have NLOS range 
measurements The standard deviation a m (m) of the smoothed curve determined in Fig 4 is 
shown in Table 2 



TABLE 2 

STANDARD DEVIATION OF MEASUREMENTS FROM SMOOTHED 
CURVE FOR FOUR NLOS MEASUREMENTS 



Base 


NLOS 


(m) 


101 


Yes 


440.2 


102 


Yes 


444.4 


103 


\es 


463.6 


104 


Yes 


450.2 



The results indicate a similar standard deviation &(m) for all four base stations 101, 
102, 103 and 104 having NLOS 

In a third example, three results were determined using = -1 18 3m y 0 = -3 7m with 
the residual analysis tracking method shown in Fig 5. In test 1, base station 104 was NLOS 
In test 2, base station 103 and base station 104 are NLOS In test 3, base station 102, base 
station 103 and base station 104 were non-line of sight. The number of times each base 
station had the largest absolute residual difference is shown in Table 3. 
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TABLE 3 

PERCENTAGE OF TIME BS HAD LARGEST RESIDUAL 



TEST 




BS101 


BS102 


BS103 


BS104 


■ 

1 


LOS 


10 


11 


18.5 






NLOS 








60 


2 


LOS 


18.5 


15 








NLOS 






26.5 


40 


3 


LOS 


12.5 










NLOS 




20 


40 5 


27 



The results indicate NLOS base stations having larger percentages of residual 
differences. 

In a fourth example, results of the method for locat.on est.mation at the present 
invention were compared with a conventional least square analys.s, a least square analysis with 
all range measurements are line of sight and a convent.onal Cramer Rao Lower Bound 
analysis The Cramer Rao Lower Bound represents a lower bound on the rms error of any 
unbiased estimator Table 4 represents the present method shown in column 2, the 
conventional least squares analysis shown in column 1, a least square analysis with all 
measurements LOS in column 3 and the conventional Cramer Rao Lower Bound analysis 
shown in column 4. The location and speed errors in each coordinates were measured in 
meters and meters/second respectively. 
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p.™ = mean error in estimating x<> 
M-yo = mean error in estimating y 0 
|ivx = mean error in estimating v x 
[i^y = mean error in estimating v y 



o xo = standard deviation of x 0 
o yo = standard deviation of y 0 
o vo = standard deviation of v 
Ovy - standard deviation of v 



TABLE 4 

COMPARISON OF ESTIMATOR PERFORMANCE 





LEAST SQUARES 
PRIOR ART METHOD 


METHOD OF PRESENT 
INVENTION 


LOS 


•JCRI.B 




297.8 


-3 98 


0 17 






32.9 


_ 28.30 


1642 


15.88 




-306 1 


-2.36 


0.54 




O vo 


55.5 


45.13 


14.15 


14 18 


Mvx 


0 18 


-0 09 


-0.005 




Ovx 


0.55 


0.49 


0.27 


0.27 


Mw 


4.49 


-01 


-0.005 






0.84 


0.64 , 


0.25 


0.25 



The results indicate that the mobile location estimation method of the present invention 
significantly reduced the estimation bias as compared to results without NLOS error 
correction. 

Fig. 12 is a comparison of the probability of detecting an NLOS range measurement 
The sampling period was 0.5 seconds. The number of samples varied between 5 and 1 50 X„ 
was 200m and y. was 100m Base station 101 and base station 104 were LOS Base station 
102 and base station 103 were NLOS The results indicate NLOS can be detected with high 
probability for a small number of samples 
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5 It is to be understood that the above-described embodiments are illustrative of only a 

few of the many possible specific embodiments which can represent applications of the 
principles of the invention Numerous and varied other arrangements can be readily devised 
accordance with these principles by those skilled in the an without departing from the spirit 
and scope of the invention 



in 
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We claim : 

1 A method for mobile station location estimation comprising the steps of: 

a obtaining range measurements between said mobile station and a base station; 

b. identifying whether said base station is line of sight with said mobile station or 
non-line of sight with said mobile station at the time at which a mobile location estimate is made; 

c. correcting non-line of sight range measurements for a base station identified as 
non-line of sight with said mobile station in step b to determine reconstructed line of sight range 
measurements; 

d. repeating steps a through c for a predetermined number of said base stations; 

and 

e. determining said mobile station location estimation from said reconstructed line 
of sight range measurements determined in step c or said range measurements determined in step 
a for an identified line of sight base station in step b. or the combination of said reconstructed line 
of sight range measurements determined in step c. and said range measurements determined in 
step a for an identified line of sight base station in step b 



2. The method of claim I wherein step b comprises the steps of: 
obtaining line of sight range measurements between said mobile station and a base 
station without noise; 

obtaining noisy line of sight range measurements between said mobile station and 
said base station. 
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predetermining a first standard deviation of the difference of said line of sight 
range measurements with said noisy line ot sight range measurements. 

smoothing said range measurements determined in step a; 

determining a second standard deviation of the difference between said smoothed 
range measurements and said noisy line oi sight range measurements, and 

discriminating between said base station being line of sight or said base station 
being non-line of sight from said first standard deviation and said second standard deviation, 
wherein said base station is determined t< be non-line of sight when said second standard 
deviation is greater than said first standard deviation and line of sight when said second standard 
deviation is on the order of said first star dard deviation 

3. The method of claim : wherein the range measurement obtained in step a is 

represented by: 

r„ (t, ) = L (/, ) + n m (/, ) + NLOS m (7, ) 
form= 1. ...,M / = 0, ... K-l, wherein 

U (/, ) is the LOS distance between a mobile station and the m ,h base station in two dimensions 
which is given by: 

L m (/,)= | v(/, ) +j *y(t,) - Xm -j *y„\- 
J = v-7 , | | is absolute value, 

*('. ),y(t, ) and {x m , y m ) are respective! the coordinates of the mobile station at time, /, , and 
those of the m' h base station; n m (t. ) represents conventional measurement noise such as additive 
white Gaussian measurement noise an NLOS. (/,) represents NLOS measurement error at time 
t, \ and M is the total number of base stations; and K is the total number of time samples 
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4. The method of claim 3 where the range measurement is smoothed by modeling: 

r m {t,) = 2>„, (")'." 

and solving for the unknown coefficients, {a m (n)) with a least squares technique. 

5. The method of claim 4 wherein the second standard deviation is represented by 



wherein 

.V-/ 

s„(t, )=£ a m (/;)/;. 

n-0 

6 The method of claim 1 wherein step b comprises the steps of: 
estimating coordinates of said mobile station from said range measurement 
obtained in step a over time; 

calculating a range measurement from said estimated coordinates, 
determining a residual from the difference of said range measurement obtained ii 
step a and said calculated range measurement, 

counting the number of times the residual is the greatest at each base station for 
each time instant, and 
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defining said base station as non-line of sight from the base station which has the 
greatest value of the number of times the greatest residual was counted 

7 The method of claim 6 wherein said estimated coordinates are represented by 
*«(',)..?«(/,) at each instance of time /, , said estimated coordinates are determined as least 

Af 

squares estimates to F = £(/■„(/, ) - L m {t t )) 2 

m.l 

where /:„(/,) = | i (/, ) - x m + j* y(t,)- j* y m \ 

8 The method of claim 1 wherein step c comprises the steps of: 

determining a value of maximum noise deviation and standard deviation from said 

range measurements obtained in step a and a predeterm.ned line of sight range measurement with 
negligible noise; 

smoothing said range measurements obtained from step a; 
graphing a curve of said smoothed range measurements; 

determining a point of maximum deviation of said range measurement below said 

curve; 

displacing said curve downwards to pass through said point of maximum 

deviation; and 

displacing said curve upwards by said value of said maximum noise deviation, 
thereby providing said reconstructed range measurement 
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9. The method of claim 1 wherein said steps a through c are repeated for at least 
two base stations and further comprising the step of determining angle arrival information. 

wherein said mobile station location is estimated from range measurements or 
reconstructed line of sight range measurements of said two base stations and said angle arrival 
information. 

10. The method of claim 1 wherein steps a through c are repeated for three base 

stations 

1 1. A system for mobile station location estimation comprising 

means for obtaining range measurements between said mobile station and a 
plurality of base stations; 

identifying means for identifying whether each of said base stations is line of sight 
with said mobile station as a line of sight base station or non-line of sight with said mobile station 
as a non-line of sight base station, 

correcting means for correcting said range measurement for each of said non-line 
of sight base stations to determine a reconstructed line of sight range measurement; and 

estimating means for determining said mobile station location estimation from said 
reconstructed line of sight range measurements or said range measurements for said line of sight 
base station, or the combination of said reconstructed line of sight range measurements, and said 
range measurements for said line of sight base stations 
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12. The system of claim 1 1 wherein said identifying means comprises: 
means for obtaining a line of sight range measurement without noise between said 
mobile station and each of said base station; 

means for obtaining a noisy line of sight range measurement between said mobile 
station and each of said base station; 

means for predetermining a first standard deviation of the difference of said line of 
sight range measurement with said noisy line of sight range measurements, 
means for smoothing said range measurements; 

means for determining a second standard deviation of the difference between said 
smoothed range measurements and said noisy line of sight range measurement; and 

means for discriminating each of said base stations as being line of sight or being 
non-line of sight from said first standard deviation and said second standard deviation, wherein 
said base station is determined to be non-line of sight when said second standard deviation is 
significantly greater than said first standard deviation and line of sight when said second standard 
deviation is on the order of said first standard deviation. 

13. The system of claim 12 wherein the range measurement is tepresented by 
r m (t, ) = L m (t, ) + n m {t, ) + NLOS m (t, ) 

for/w= 1, ...,M / = 0, ... K-l, wherein 

U (/, ) is the LOS distance between a mobile station and the m* base station in two d.mensions 
which is given by: 

L m (l, ) = \x(t, ) +j . >•(,, ) - Xm - J , y m | . 
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is absolute value. 

x (t>),Xt> ) and (x m ,y m ) are respectively the coordinates of the mobile station at time, /, , and 
those of the m' h base station; n m (t ; ) represents conventional measurement noise such as additive 
white Gaussian measurement noise and NLOS m (/,) represents NLOS measurement error at time 
/,.; and Mis the total number of base stations; and K is the total number of time sample 



les. 



14. The system of claim 13 where the range measurement is smoothed by 

modeling: 



r„(/.) = £> m (/;)/; 

n = 0 

and solving for the unknown coefficients, {a m («)} ^ with a least squares techniq 



ue. 



15 The system of claim 14 wherein the second standard deviation is represented 

by 



//*'-' 



wherein 



N-l 



16. The system of claim 12 wherein said identifying means comprises: 
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means for estimating coordinates of said mobile station from said range 
measurements from a plurality of base stations received over time; 

means for calculating a calculated range measurement from said estimated 

coordinates; 

means for determining a residual from the difference of said range measurements 
and said calculated range measurement, 

means for counting he number of times the residual is the greatest at each base 
station for each time instant; and 

means for defining aid base station as non-line of sight from the base station 
which has the greatest value of the number of times the greatest residual was counted 

17. The system ot claim 14 wherein said estimated coordinates are represented by 
x is (',).>'«('.) at each instance ,f time /, , said estimated coordinates are determined as least 



M 



squares estimates to F = £ (r„ ( /, ) - L m ( t, )) 

where L m (t,)= I i(t,)-x m +;^(f,)- y. y m \ 

18. The system <-f claim 12 wherein said estimating means comprises: 

means for determining a value of maximum noise deviation and standard deviation 

for each of said range measurer >ents and a predetermined line of sight range measurement with 

negligible noise; 

means for smoothing said range measurements; 



26 



WO 98/14018 



PCT7US97/17786 



means for graphing a curve of said smoothed range measurements; 
means for determining a point of maximum deviation of said range measurements 
below said curve; 

means for displacing said curve downwards to pass through said point of 
maximum deviation; and 

means for displacing said curve upwards by said value of said maximum noise 
deviation, thereby providing said reconstructed range measurement 

19. The system of claim 12 further comprises means for obtaining angie arrival 
information wherein said mobile station location is estimated from range measurements or 
reconstructed line of sight range measurement of said base stations and said angle arrival 
information 



27 



WO 98/14018 



PCT/US97/17786 




1/13 



WO98/14018 



PCT7US97/ 17786 




WO 98/14018 



PCT/US97/ 17786 



20 



21 

\ 



Obtain range measurements 



22 



IDENTIFY IF BASE STATION IS 
LOS ORNLOS 



24 

\ 



26 

\ 



NLOS 


FOR 

\ 


CORRECT NLOS RANGE MEASUREMENT 
TO RECONSTRUCT AS AN LOS 




RECONSTRUCTED 
LOS 

I 


DETERMINE MOBILE LOCATION 
ESTIMATION 



1 



1...M 



LOS 



FIG. 2 



3/13 



WO 98/14018 



PCT/US97/17786 




4/13 



WO 98/14018 



PCT/US97/17786 



30 

^ Obtain LOS range measurement between base 
station and mobile station 



i 

Obtain noisy LOS range measurement 
between base station and mobile station 



32 

\ 



I 

Predetermine standard deviation and 
maximum deviation of LOS range 
measurement with noisy LOS range 
measurement 



22 



33 



34 

\ 



Smooth range measurement obtained from 
block 21 



1 



Determine standard deviation of the difference 
between smoothed range measurement and noisy 
range measurement (i.e., standard deviation of 
residual) 



I 



^ From value of standard deviation from block 32 
and standard deviation from block 34 
discriminate between LOS or NLOS range 
measurements 



FIG. 4 



5/13 



PCT/US97/17786 



Obtain range measurements between mobile 
and base station 



i 

Use all range measurements to estimate 
mobile station coordinates over time 



i 

Calculate a range from the estimated coordinates 
and known base station coordinates 



Determine residual as difference between 
measured range and calculated range using 
estimated coordinates 



i 

Count the number of times the residual is the 
greatest at each base station for each time 
instant, /, 



From value of residual differences determine 
which base station has a significantly larger 
number of greatest absolute residual differences 
than the number of greatest absolute residual 
differences from other base stations and declare 
it NLOS 



FIG. 5 

6/13 



WO 98/14018 



PCT/US97/17786 



Smooth range measurements obtained from 
block 21 



1 

Graph curve of smoothed range 
measurements 



^ Determine point of maximum deviation of the 
noisy range measurement containing the 
NLOS error below the smoothed curve 



i 

Displace smoothed curve downwards to the 
point of maximum deviation so that the 
smoothed curve passes through the point of 
maximum deviation 



Displace curve upwards by value of noise 
deviation, cc m for reconstructing measurement 
ofLOS 



FIG. 6 



7/13 



WO 98/14018 



PCT/US97/17786 



8S00 



6000 



7500 



7000 



E 6500 - 



o> 
c 



£ 6000 



5500 



5000 



4500 



Corrected Range 

True Range 

- - Smoothed Range 



4000 



10 




20 



30 



40 50 60 
Time (s) 



70 



i 

80 



90 100 



AD 



8/13 



WO 98/14018 



PCT/US97/17786 




9/13 



WO 98/1 4018 PCTYUS97/1 7786 




10/13 



WO 98/14018 



PCT/US97/17786 



f [isr ID ft 

2-D Tracking Error without NLOS Identification and Correction 
-150r 1 1 1 1 , , 



-200 - 
-250 - 
-300- 




-350' 1 1 1 " 1 1 1 

80 100 120 140 180 180 200 220 240 

x-coordinate (m) 

2-D Tracking Error after NLOS Identification and Correction 

30 1 1 r— — i i r 1 1 1 




-20 ' 1 ! 1 1 1 i — i 

-50 -40 -30 -20 -10 0 10 20 30 

x-coordinate (m) 



11/13 



WO 98/14018 



PCT/US97/ 17786 



f\*r II A 

2-D Tracking Error without NLOS Identification and Correction 




200 



-40 ' 1 " ! — -J L 

^60 -40 -20 0 20 40 60 

x-coordinate (m) 

-P l(r IIS 



12/13 



WO 98/14018 PCTYUS97/I7786 



0.9 




Number of Samples Per Trial 



Pi 0 i >- 



13/13 



INTERNATIONAL SEARCH REPORT 



International application No. 

PCT/US97/17786 



A. CLASSIFICATION OF SUBJECT MATTER 
1PC(6) :H04Q 07/20 
US CL :455/456 

According to International Patent Classification (IPC) or to both national classification and IPC 



FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

U.S. : 455/456, 403. 422. 435. 457. 504. 506. 550; 701/214. 215; 342/450. 357. 457 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched 
NONE 



Electronic dau base consulted dunng the international search (name of data base and. where practicable, search terms used) 
APS. IEEE 

search terms: LOS. NLOS. line of sight, multipath fading, standard deviation, curve fit, smoothing. 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 4 



Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



Y, P 
Y 



US, 5,629,707 A (HEUVEL et al) 13 May 1997, see Col. 6;55 - 
Col.l0;65. 

US 5,257,405 A (REITBERGER) 26 October 1993, see Col. 9;55 - 
Col. 10;56. 

US 5,327,144 A (STILIP et al) 05 July 1994, see CoL 14;15-45, 
CoL 16;5 - Col. 17;65. 

US 5,390,124 A (KYRTSOS) 14 February 1995, see whole 
document. 

BORJESON et al, Outdoor Microcell Measurements at 1700 MHz, 
IEEE Vehicular Technology Conference, 1992, pp. 927-931 



1-2, 6, 8-12, 16, 
18-19 

1-2, 6, 8-12, 16, 
18-19 

1-2, 6, 8-12, 16, 
18-19 

1-19 



1-19 



[ 1 Farther documents are listed in the coofaauatioo of Box C. 



patent family annex. 



loWof paniwsf ftoranaa 



1 fiUnt data or priority 



l vtbofc mmy daw 4oukm am priority i1mi(i) or vkte* i 



«alor ■■■■ tbi 

MB«BI M tokao i 



to oroKt an arrmotrw step 



i (as apactAad) 



poor tot 
afcamad 



doauaaa* of parttouter ratevanoa; lha oUn*d brv«ot> 
aaaaaaarad to awoHa mm a»aa*nra atop vbao ttto 



fdtoa data but later than 



i to a par* 
aibar of* tba aaaia patent family 



Dam of me actual completion of tac international 
09 DECEMBER 1997 



Date of mailing of tbe international search report 




1 0 FEB 1998 



Name and mailing address of the ISA/US 
Coasminioeer of Patents and Tn 
BoxfCT 

WsUuDgtoo, D.C 30231 
Facsimile No (703) 305-3230 



,1AM TROST 
Telephone No. (703) 30S-S3 1 8 



Form PCT/I&A/2I0 (second sneetXJaly 1992)* 



CORRECTED 
VERSION* 



WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau 




PCT 

INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) International Patent Classification 6 
H04Q 7/20 



Al 



(11) International Publication Number: 
(43) International Publication Date: 



WO 98/14018 

2 April 1998 (02.04.98) 



(21) International Application Number: PCT/US97/I7786 

(22) International Filing Date: 29 September 1997 (29.09.97) 



(30) Priority Data: 

60/027.453 



27 September 1996 (27.09.96) US 



(71) Applicant: RUTGERS UNIVERSITY [US/US]; Office of 

Corporate Liaison & Technology Transfer, Administrative 
Services Building, Annex II, Piscataway, NJ 08854-1 179 
(US). 

(72) Inventors: WYLIE, Marilynn, P.; 1621 North Robinson Road, 

Philadelphia, PA 19151 (US). HOLTZMAN, Jack, M.; 
Apartment 2A, 10 West Court, Westfiled, NJ 07090 (US). 

(74) Agent: McKAY, Diane, Dunn; Mathews, Collins. Shepherd & 
Gould, Suite 306, 100 Thanet Circle, Princeton, NJ 08540 
(US). 



(81) Designated States: AL, AM, AT, AU, AZ. BA, BB, BG. BR, 
BY, CA, CH, CN, CU, CZ, DE, DK, EE, ES, FI, GB, GE, 
GH, HU, ID, IL, IS, JP. KE, KG, KP, KR, KZ, LC, LK, 
LR, LS. LT, LU, LV, MD, MG, MK, MN, MW, MX, NO. 
NZ. PL, PT, RO, RU, SD, SE, SG, SI, SK, SL, TJ, TM, 
TR, TT, UA, UG, UZ, VN, YU, ZW, ARIPO patent (GH, 
KE, LS, MW, SD, SZ, UG, ZW), European patent (AT, BE, 
CH, DE, DK, ES, FI, FR. GB, GR, IE, IT, LU, MC, NL, 
PT. SE), OAPI patent (BF, BJ, CF, CG, CI, CM, GA. GN, 
ML, MR, NE, SN, TD, TG). 



Published 

With international search report. 

Before the expiration of the time limit for amending the 
claims and to be republished in the event of the receipt of 
amendments. 



(54) Title: METHOD AND SYSTEM FOR MOBILE LOCATION ESTIMATION 



(57) Abstract 

The present invention relates to a method 
and system for mobile location estimation in which 
range measurements between the mobile station 
and a base station are discriminated as being ei- 
ther from a base station which is line of sight 
with the mobile station or a base station that is 
non-line of sight with the mobile station (22). 
Non— line of sight error present in range measure- 
ments from base stations that are non-line of sight 
with the mobile station is corrected (24). A recon- 
structed line of sight range measurement is estab- 
lished from the correction of the non-line of sight 
error. The reconstructed line of sight range mea- 
surements can be used with range measurements 
of base stations determined to be line of sight for 
accurately locating the mobile station (26). 



OBTAIN RANGE MEASUREMENTS 



■21 



IDENTIFY IF BASE STATION 
IS LOS OR NLOS 



FOR 1....M 



NLOS \3 



22 



24 



CORRECT NLOS RANGE MEASUREMENT 
TO RECONSTRUCT AS AN LOS 



26 



RECONSTRUCTED 
LOS 



DETERMINE MOBILE 
LOCATION ESTIMATION 



LOS 



FOR THE PURPOSES OF INFORMATION ONLY 



Codes used to identify States party to the PCT on the front pages of pamphlets publishing international applications under the PCT. 



AL 


Albania 


ES 


Spain 


IS 


Lesotho 


SI 


Slovenia 


AM 


Armenia 


FI 


Finland 


LT 


Lithuania 


SK 


Slovakia 


AT 


Austria 


FR 


France 


LU 


Luxembourg 


SN 


Senegal 


AU 


Australia 


GA 


Gabon 


LV 


Latvia 


sz 


Swaziland 


AZ 


Azerbaijan 


GB 


United Kingdom 


MC 


Monaco 


TD 


Chad 


BA 


Bosnia and Herzegovina 


GE 


Georgia 


MD 


Republic of Moldova 


TG 


Togo 


BB 


Barbados 


GEI 


Ghana 


MC 


Madagascar 


TJ 


Tajikistan 


BE 


Belgium 


GN 


Guinea 


MK 


The former Yugoslav 


TM 


Turkmenistan 


BF 


Burkina Faso 


GR 


Greece 




Republic of Macedonia 


TR 


Turkey 


BG 


Bulgaria 


HU 


Hungary 


ML 


Mali 


TT 


Trinidad and Tobago 


BJ 


Benin 


IE 


Ireland 


MN 


Mongolia 


UA 


Ukraine 


BR 


Brazil 


IL 


Israel 


MR 


Mauritania 


UG 


Uganda 


BY 


Belarus 


IS 


Iceland 


MW 


Malawi 


US 


United States of America 


CA 


Canada 


IT 


Italy 


MX 


Mexico 


UZ 


Uzbekistan 


CF 


Central African Republic 


JP 


Japan 


NE 


Niger 


VN 


Viet Nam 


CG 


Congo 


KE 


Kenya 


NL 


Netherlands 


YU 


Yugoslavia 


CH 


Switzerland 


KG 


Kyrgyzstan 


NO 


Norway 


ZW 


Zimbabwe 


CI 


Cdte d'lvoire 


KP 


Democratic People's 


NZ 


New Zealand 






CM 


Cameroon 




Republic of Korea 


PL 


Poland 






CN 


China 


KR 


Republic of Korea 


PT 


Portugal 






cu 


Cuba 


KZ 


Kazakstan 


RO 


Romania 






cz 


Czech Republic 


LC 


Saint Lucia 


RU 


Russian Federation 






DE 


Germany 


LI 


Liechtenstein 


SD 


Sudan 






DK 


Denmark 


LK 


Sri Lanka 


SE 


Sweden 






EE 


Estonia 


LR 


Liberia 


SG 


Singapore 







WO 98/14018 



PCT/US97/17786 



METHOD AND SYSTEM FOR MOBILE LOCATION ESTIMATION 



This application claims the benefit of U.S. Provisional Application No. 60/027,453 
entitled Non-Line Of Sight Problem in Mobile Location Estimation filed by Applicants on 
September 27, 1996 hereby incorporated by reference into this application. 

Background of the Invention 

1 Field of the Invention 

The present invention relates to a method and system for mobile station location 
estimation in which base stations that are in line of sight of the mobile station and base stations 
that are not in the line of sight of the base station can be determined. Errors in base station 
signals generated from determined non-line of sight base stations are reduced for providing 
improved mobile station location estimation. 

2 Description of the Related Art 

Mobile location estimation determines a geographical estimate of the location of a 
mobile station. Mobile location estimation is useful in management of fleets of mobile 
stations, location dependent information services, location dependent billing services and 
Emergency 911 location of a mobile station Enhanced 911 is designed to automatically 
forward the number of a caller to a public safety answering point (PSAP) In implementing 
enhanced 911 in a wireless network, wireless service providers provide two d.mensional 
location of the vehicle to the public safety answering point (PSAP). The Federal 
Communications Commission (FCC) has regulated by the year 2001 that wireless service 
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5 providers have the capability of locating callers in two dimensions within 125 meters 67% of 
the time. 

One conventional method for locating a mobile station in two dimensions would use 
the measurement of the line of sight distance between the mobile station and at least three base 
stations. U.S. Patent No. 5,365,516 describes a method for determining the location of a 

10 transponder unit in which a radio signal is sent by the mobile station. The arrival time of the 
radio signal is measured at each of three base stations. Each distance measurement between 
the mobile station and one of the base stations can be used to generate a circle which is 
centered at the measuring base station. The circle has a radius which is equal to the distance 
between the mobile station and the base station. Accordingly, three circles are generated, one 

15 for each of the base stations. In the absence of any measurement error of the distance between 
the base stat.ons and the mobile station, the intersection of the three circles unambiguously 
determines the location of the mobile station. This method has the drawback that the distance 
measurements can be corrupted by noise resulting in errors in determining the location of the 
mobile station. 

20 A conventional solution for providing more accurate position estimates is to reduce the 

error due to noise with a least squares analysis. Accordingly, the least squares analysis 
provides a more accurate position estimate. This solution has the limitation of not accounting 
for the possibility of a lack of a direct path between the base station and the mobile station. 
For example, in an urban environment, a building or buildings may be in the path between the 

25 mobile station and the base station. A propagating signal between the mobile station and the 
base station can be reflected and defracted by the object in the path of the mobile station to the 
base station resulting in the signal traveling excess path lengths. The excess path lengths can 
be on the order of a hundred meters 
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5 The lack of direct path between the mobile station and the' base station can be defined 

as a non-line of sight (NLOS). The importance of detecting and reducing the NLOS 
measurements between a mobile station and a base station is recognized in M.I. Silventoinen, 
et al., "Mobile Station Locating in GSM', IEEE Wireless Communication System S^nposjum 
Long Island NY, November 1995 and J.L. CafFrey et al., "Radio Location in Urban CDMA 
10 Microcells", Proceedings of the Personal Indoor and Mo bile Radio Fnvirnnm^t 1995 

U.S. Patent No. 5,365,516 ('516 patent) describes an embodiment of a transceiver 
locating system operating in an environment susceptible to multipath interference. The system 
includes a transponder which is operable within a prescribed coverage area to transmit a burst 
of data symbols in a coded carrier pulse. Each base station includes a receiver for detecting 
15 and responding to the data symbol at a given time, interrupting the data symbol and rejecting 
echoes resulting from multipath interference. A comparison circuit responds to the receiver 
for comparing respectively identified given times and decorrelating the time difference to 
improve data quality. Although the '516 patent addresses multipath interference, it does not 
attempt to detect base stations for reducing multipath NLOS with mobile stations 
20 It is desirable to provide a method and system for providing improved mobile location 

estimation which is robust to NLOS error. 

Summary of the Invention 
Briefly described, the present invention relates to a method and a system for mobile 
location estimation in which base stations are identified to be either line of sight (LOS) or non- 
25 line of sight (NLOS) with a mobile station. A range measurement is determined as the distance 
between the base station and the mobile station. NLOS ranging error is corrected for base 
stations identified to be NLOS with the mobile station by reconstructing the LOS 
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5 measurement. From the range measurements of base stations identified as LOS and the 
reconstructed LOS range measurements the location of the mobile station is estimated. 

The base station can be identified as being NLOS by comparing the standard deviation 
of standard measurement noise from the environment to the standard deviation of a smoothed 
range measurement obtained from the range measurements between the base station and 

10 mobile station. The smoothed range measurement can be obtained using an N* order 
polynomial fit. It has been found that when the standard deviation of the smoothed range 
measurement is on the order of the standard deviation of the standard measurement noise, the 
base station corresponds to an LOS environment and when the standard deviation of the 
smoothed range measurement is greater than the standard deviation due to standard 

1 5 measurement noise, the base station corresponds to an NLOS environment. Alternatively, the 
residuals from a least squares analysis can be used to determine the presence of NLOS range 
measurements. 

NLOS error can be corrected when the standard measurement noise dominates the 
NLOS error and there is predetermined identification of the approximate support of the 

20 standard measurement noise over the real axis. A reconstructed LOS range measurement can 
be determined by graphing a curve of the smoothed range measurements The point of 
maximum deviation of the smoothed range measurement below the curve is determined. The 
curve is displaced downwards to pass through the point of maximum deviation. Thereafter, 
the curve is displaced upwards by the value of the maximum standard measurement noise 

25 deviation from an LOS measurement with negligible noise, thereby providing a reconstructed 
range measurement. 

The mobile location estimation can be determined using at least three range 
measurements between LOS base stations and the mobile station or reconstructed LOS range 
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measurements in a multilateration analysis. In this analysis, a circle is generated from each 
range measurement. The circle is centered at the base station and the range measurement is 
the radius of the circle. The estimated intersection of the three circles determines the location 
of the mobile station. Alternatively, two range measurements and information directed to the 
position angle of the mobile station can be used for estimating the location of the mobile 
station. 

The present invention has the advantages of accurately determining the location of a 
mobile station by reducing NLOS error. In addition, the present invention can provide 
confidence ir; an LOS environment that all base stations are LOS with the mobile station 
Results indicate that position range bias due to NLOS error can be reduced several orders of 
magnitude with the method of the present invention. 

The present invention will be more fully described by reference to the following 
drawings. 

Brief Description of the Drawing s 

Fig. 1A is a schematic diagram of an environment in which there is an unobstructed 
line of sight radio signal path between a mobile station and a base station 

Fig. IB is a schematic diagram of an environment in which there is a non-line sight 
radio signal path between a mobile station and a base station. 

Fig. 2 is a flow diagram of the system and method for mobile location estimation in 
accordance with the teachings of the present invention. • 

Fig. 3 is a schematic diagram of distance measurements of a reconstructed line of sight 
base station and determined line of sight base stations. 
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5 Fig. 4 is a flow diagram of a method for identifying non-line of sight base stations of 

the present invention. 

Fig. 5 is a flow diagram of an alternate method for identifying non-line of sight base 
stations. 

Fig. 6 is a flow diagram of a method for reconstructing a line of sight base station for 
10 non-line of sight measurements. 

Fig. 7 is a graph of a comparison of NLOS measurements and reconstructed LOS 
measurements. 

Fig. 8 is a schematic diagram of a system for implementing the method of the present 
invention 

15 F 'S 9 is a schematic diagram of positioning of base stations used in examples of 

performance of the method of the present invention. 

Fig. 10A is a graph of two dimensional tracking without non line of sight error 
detection and correction. 

Fig. 1 OB is a graph of two dimensional tracking with non line of sight error detection 
20 and correction. 

~- Fig. 1 1 A is a graph of two dimensional tracking without non line of sight error 
detection and correction. 

Fig. 1 IB is a graph of two dimensional tracking with non line of sight error detection 
and correction. 

25 Fi 8 !2 is a graph of the fraction of time a base station was declared NLOS using the 

residual rank analysis method 
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5 Detailed Description of the Present Invention 

During the course of this description like numbers will be used to identify like elements 
according to the different figures which illustrate the invention. 

Fig. 1A illustrates a schematic diagram of a line of sight (LOS) path 10 between a base 
station 12 and mobile station 14. Signal 13 can be transmitted from base station 12 to mobile 
10 station 14 and returned from mobile station 14 to base station 12 Fig. IB illustrates a 
schematic non-line of sight (NLOS) path 1 1 between base station 12 and mobile station 14. 
Building 15 is positioned between base station 12 and mobile station 14 resulting in reflection 
of signal 16. For example, signal 13 and signal 16 can be a radio signal 

A range measurement for measuring the distance between base station 1 2 and mobile 
15 station 14 can be measured as the time it takes a signal sent between base station 12 and 
mobile station 14: 

r = cT (1) 
in which the mobile station to base station range measurement is represented by r, c represents 
the speed of light which is the same speed as the propagation of radio waves and T represents 
20 the one-way travel time of the signal. A range measurement of the distance between mobile 
station 14 and base station 12 in Fig. 1 A and IB can be determined using equation (1) based 
on travel time of signal 13 and signal 16, respectively, between base station 12 and mobile 
station 14. The value of r generated from signal 16 is greater than the value of r generated for 
signal 13. 

25 F 'g- 2 is a flow diagram of the system and method of the present invention for mobile 

location estimation 20. In block 21, a range measurement is obtained between mobile station 
14 and base station 12 using equation (1). In block 22, base station 12 is identified as being in 
line of sight (LOS) or non-line of sight (NLOS) with mobile station 14. Block 22 is repeated 
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for a plurality of base stations 12 positioned at different locations from mobile station 14. If 
base station 12 is identified to be LOS in block 22, the range measurement obtained from 
block 21 is forwarded to block 26. If base station 12 is identified to be NLOS in block 22, 
block 24 is implemented for reducing the error of the range measurement between base station 
12 and mobile station 14, thereby rendering the range measurement between base station 12 
and mobile station 14 as a reconstructed LOS base station 13, as shown in Fig. 3. 

In Fig. 3, base station 12 labeled BS1 has a range measurement labeled LOS RANGE 
1 determined to be LOS. Base station 12 labeled BS2 has a range measurement labeled LOS 
RANGE 2 determined to be LOS. Base station 13 labeled BS3 has a range measurement 
labeled NLOS RANGE 3 determined to be NLOS. A range measurement for the 
reconstructed LOS base station labeled RECONSTRUCTED RANGE 3 is forwarded to block 
26. Range measurements from determined LOS base stations from block 22 labeled LOS 
RANGE 1, and LOS RANGE 2 are also forwarded to block 26. From the range 
measurements of the reconstructed LOS base stations or the determined LOS base stations, or 
a combination of range measurements of the reconstructed LOS the mobile location estimation 
can be identified using a conventional multilateration technique, such as described in U.S. 
Patent No. 5,365,516, hereby incorporated by reference into this application Alternatively, 
the mobile location estimation can be determined from time difference of arrival time 
measurements as the difference of propagation delays between the mobile station 14 and pairs 
of base stations 12. In this case, the position estimate is at the intersection of hyperbolas. The 
number of base stations can be reduced below three if there is also angle of arrival information. 
These methods are described in IS. Rappaport et al., " Position Location Using Wireless 
Communication On Highways Of the Future", IEEE Commun ications Marine October 
1996. 
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One method for identifying if a base station is LOS or NLOS in block 22 is illustrated 
in Fig. 4. In this method, the time history of range measurements between base station 12 and 
mobile station 14 is combined with predetermined standard deviation from conventional 
measurement noise in a radio signal environment. 

The arrival time of signals sent from base station 12 to a mobile station 14 and 
transponded back to the base station 12 can be converted to a range measurement, in block 
30. The range measurement at the m* base station at time t k can be represented as: 
r m (/, ) = L m (/, ) + n m (/, ) + NLOS m (/, ) (2) 

for m = 1, M / = 0, ... K-l, wherein 

L m (/, ) is the LOS distance between a mobile station and the base station in two 
dimensions which is given by: 

L m (t,)= \x(t<)+j *y(i,) -x m -j *y m \\ (3) 
)> y{*i ) and (x m , y m ) are respectively the coordinates of the mobile station at time, l, , and 
those of the m th base station; n m (t t ) represents conventional measurement noise such as 
additive white Gaussian measurement noise and NLOS™ (/,) represents NLOS measurement 
error at time t,.\ and M is the total number of base stations; and A' is the total number of time 
samples. 

In block 30, an LOS range measurement with negligible noise is obtained for base 
station 12 in LOS with mobile station 14. The LOS range measurement can be obtained by 
physically measuring a range between base station 12 and mobile station 14 or can be obtained 
as a range measurement determined by equation (1) in a negligible noise environment. In 
block 31, a noisy range measurement is determined as a range measurement which is LOS 
with a base station taken in a noisy environment In block 32, the standard deviation of the 
noisy range measurement from the LOS measurement without noise is determined. Blocks 30, 



31 and 32 can be predetermined before identifying base station 12 as either LOS or NLOS in 
block 22. The standard deviation due to noise n m (() can be represented by o m . 

In block 32, the range measurement obtained from block 21 is smoothed by modeling 

r m (/ # ) = Xa m (/,)/; (4) 

and solving for the unknown coefficients, {a m (n)}^ with a least squares technique. The 
smoothed range measurement can be represented as: 

s m (t,)=^l a m (n)r:. (5) 

n-0 

In block 34, the standard deviation of the smoothed range measurement from a noisy 
range measurement (i.e., the residual) is determined. The standard deviation of the residual 
from block 34 can be represented as a m since a; = Efn 2 m (/)}. The smoothed range 
measurements along with the noisy range measurement can be used to determine standard 
deviation a m with the formulation of: 

° m = ^^(sM-rM) 2 (6) 

From the value of the standard deviation, 6 m and the standard deviation a m , the range 
measurement can be determined as either the result of base station 12 being LOS or NLOS, in 
block 36. When the range measurement has NLOS error, the value of the standard deviation 
a m is significantly larger than the value of the standard deviation a w . Accordingly, range 
measurement for base station 12 that is NLOS with mobile station 14 is determined when the 
& m is greater than the standard deviation o„, . A range measurement of base station 12 that is 
LOS with mobile station 14 is determined when the standard deviation a m is on the order of 
the standard deviation o m . 
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Alternatively, a residual analysis ranking method can be used to identify a range 
measurement as being from a base station 13 NLOS with mobile station 14. Range 
measurements between mobile station 12 and base station 14 which have been obtained in 
block 21 are inputted to block 41. At each instance of time /„ estimated coordinates 
*Ls( / ,-)>.Pis(0 of mobile station 14 are determined as least squares estimates in block 41. 
The estimated coordinates 

XlsO. X.Pls(', )) are selected to minimize the formulation: 



f, :=!>,,(',)- 4(0)*' 



(7) 



where L m {t, ) = \x(t,)-x m + j* y{t,)- j* y m | . 

In block 41, a calculated range measurement is determined from the estimated 
coordinates. In block 42, a residual difference of the range measurement between mobile 
station 12 and base station 14 with the calculated range measurement is determined. The 
residual difference can be represented as: 

e m (t, ) = r m (t,)- L m (t, ) (8) 
In block 44, the number of times the residual difference of a range measurement to a 
base station 12 has the largest value in comparison to the residual difference determined for 
range measurements at other base stations is counted for each time instant t,. It has been 
found that base stations having a range measurement between a base station NLOS with a 
mobile station have a significantly larger number of greatest absolute residual differences than 
the number of greatest absolute residual differences from other base stations From the value 
of the counted number of residual differences, base station 14 can be defined as a base station 
12 that is a LOS or a base station 12 that is NLOS with mobile station 14. 
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Fig. 6 illustrates a method for correction of range measurements between a base 
station 12 that has been determined NLOS with mobile station 14 to reconstruct a LOS range 
measurement. Data related to the range measurements from block 21 are smoothed using an 
N* order polynomial fit described in block 32. The smooth range measurements are inputted 
to block 52. The maximum deviation below the smoothed curve due to NLOS error in 
determined in block 56. It has been found that NLOS error is a non-negative random variable 
which can be approximately represented in a real axis as follows. 

0 <_ NLOSm (l, ) < Pm 

in which Pm is the maximum value of NLOS error. The standard measurement no.se, n m ( tj ) 
can be represented as a zero-mean random variable which can be approximately represented in 
a real axis as follows: - a m ^n m ( ti ) < a m , so that in a range measurement in which there is also 
an NLOS error, the total noise component can be approximated represented over the real axis 
as follows: 

-a m s n m (/) + NLOS m (/) <0 m -a m . 

It has been found that the point of maximum deviation of the measured range below 
the smoothed curve is about a m below the LOS function represented as L m (t,). In block 58, 
the smoothed curve is displaced mathematically downward to the point of maximum deviation. 
The smoothed curve is displaced mathematically upward by a value of the noise deviation a m 
in block 60 to provide a reconstructed curve representing a reconstructed. LOS base station. 

Fig. 7 represents a graph of a comparison of simulated range measurements. Curve 90 
represents the true time range measurement between a base station 12 which is LOS with a 
mobile station 14. Curve 91 represents determined range measurements having NLOS error. 
Curve 92 represents a smoothed range measurement of block station 12 and mobile station 14 
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determined from block 30 of Fig. 4 Curve 93 represents base station 12 which is 
reconstructed LOS with the mobile station 14 from block 60 of Fig. 6. 

Fig. 8 is a schematic diagram of a system 80 for implementing the method for mobile 
location estimation. System 80 includes base station server 81 Base station server 81 can be a 
computer located at base station 12 or networked thereto Base station server 81 
communicates with base station 12 for requesting and receiving data related to range 
measurements of mobile station 14 and base station 12. Base station server 81 also collects 
information on range measurements between mobile station 14 and each of base stations 81A- 
81N. The information is reported to base station server 81 either by mobile station 14 or base 
station servers 81 A-81N. The functions of modules shown in Figs. 4-6 which are coded with a 
standard programming language, such as C" programming language. The coded modules can 
be executed by base station server 8 1 . 

Results for examples of mobile location estimates with system 80 are shown in Tables I 
- IV and Fig. 9 through Fig. 12. In all of the examples, the vehicle's position in the x-y plane 
at any is given by: 

x(t) = x a + vj 

>(') = y a + \\t 
x(t) represents the x - coordinate in x -y plane at time instant, t, 
y(t) represents the y - coordinate in x - y plane at time instant, t, 
x 0 represents the initial x - coordinate, 
y a represents the initial y - coordinate, 
v x represents the speed in x - direction, 
v y represents the speed in y - direction. 
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The sampling period was chosen to be 0.5s and 200 samples' were taken. The velocity 
remained constant at v s = 9.7 m/s and v y = 16.8 m/s. Base stations 12 were assigned to have 
NLOS or LOS range measurements. Standard deviation of the standard measurement noise 
was represented as a o m was 150m and B m was chosen as 1300m. In each example three base 
stations 101, 102, 103 were used uniformly spaced around a circle of 5 kilometers and a fourth 
base station 104 was located at the center of the circle, as shown in Fig. 9. 

In a first example, base station 101 and base station 102 provide NLOS range 
measurements and base stations 103 and base station 104 provide LOS range measurements. 
The standard deviation & m (m) of the smoothed curve determined in Fig. 4 is shown in Table 
1. 

TABLE 1 

STANDARD DEVIATION OF MEASUREMENTS FROM SMOOTHED 
CURVE FOR 2 NLOS MEASUREMENTS 



Base 


NLOS 


o"™ (m) 


101 


Yes 


467.3 


102 


Yes 


447.6 


, 103 


No 


163.1 


104 


No 


142.1 



The results indicate base stations 101 and 102 have NLOS range measurements with a 
significantly larger standard deviation than base station 103 and base station 104 having a LOS 
range measurement 

Fig. 10A shows two-dimensional tracking error without NLOS identification and 
correction Fig. 10B shows two dimensional tracking error after the method of mobile location 
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estimation of the present invention is performed. The results indicate improvement of 
estimated vehicle trajectory after NLOS identification and correction. 

In a second example, base stations 101, 102, 103 and 104 have NLOS range 
measurements. The standard deviation a m (m) of the smoothed cuwe determined in Fig. 4 is 
shown in Table 2. 



TABLE 2 

STANDARD DEVIATION OF MEASUREMENTS FROM SMOOTHED 
CURVE FOR FOUR NLOS MEASUREMENTS 



Base 


NLOS 


° n, (m) 


101 


Yes 


440.2 


102 


Yes 


444.4 


103 


Yes 


463.6 


104 


Yes 


450.2 | 



The results indicate a similar standard deviation &(m) for all four base stations 101, 
102, 103 and 104 having NLOS. 

In a third example, three results were determined using x<, = -1 18.3m y 0 = -3.7m with 
the residual analysis tracking method shown in Fig 5. In test 1, base station 104 was NLOS. 
In test 2, base station 103 and base station 104 are NLOS. In test 3, base station 102, base 
station 103 and base station 104 were non-line of sight. The number of times each base 
station had the largest absolute residual difference is shown in Table 3. 
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TABLE 3 

PERCENTAGE OF TIME BS HAD LARGEST RESIDUAL 



TEST 




R*\10I 

OJIUl 


dc i no 


BS103 


BS104 


1 


LOS 


10 


11. 


18.5 






NLOS 








60 


2 


LOS 


18.5 


15 








NLOS 






26.5 


40 


3 


LOS 


12.5 










NLOS 




20 


40.5 


27 



The results indicate NLOS base stations having larger percentages of residual 
differences. 

In a fourth example, results of the method for location estimation at the present 
invention were compared with a conventional least square analysis, a least square analysis with 
all range measurements are line of sight and a conventional Cramer Rao Lower Bound 
analysis. The Cramer Rao Lower Bound represents a lower bound on the rms error of any 
unbiased estimator. Table 4 represents the present method shown in column 2, the 
conventional least squares analysis shown in column 1, a least square analysis with all 
measurements LOS in column 3 and the conventional Cramer Rao Lower Bound analysis 
shown in column 4. The location and speed errors in each coordinates were measured 
meters and meters/second respectively 



in 
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|i xo = mean error in estimating x<, 
|iyo = mean error in estimating y Q 
^vx = mean error in estimating v x 
^vy = mean error in estimating v y 



PCT/US97/17786 

o xo = standard deviation of x 0 
o yo = standard deviation of y o 
o vo = standard deviation of v x 
Ovy = standard deviation of v v 



TABLE 4 

COMPARISON OF ESTIMATOR PERFORMANCE 





LEAST SQUARES 
PRIOR ART METHOD 


METHOD OF PRESENT 
INVENTION 


LOS 


ylCRLB 




297.8 


-3.98 


0.17 




Oxo 


32.9 


_ 2830 


16.42 


15.88 


Hvo 


_ -306. 1 


-2.36 


0.54 




Oxo 


55.5 


45.13 


14.15 


14.18 


Hvx 


0.18 


-0.09 


-0.005 




Ovx 


0.55 


049 


0.27 


0.27 


Uw 


4.49 


-.01 


-0.005 




Cw 


0.84 


0.64 


0.25 


0.25 



The results indicate that the mobile location estimation method of the present invention 
significantly reduced the estimation bias as compared to results without NLOS error 
correction 

Fig. 12 is a comparison of the probability of detecting an NLOS range measurement. 
The sampling period was 0.5 seconds. The number of samples varied between 5 and 150 X* 
was 200m and y c was 100m. Base station 101 and base station 104 were LOS. Base station 
102 and base station 103 were NLOS. The results indicate NLOS can be detected with high 
probability for a small number of samples. 
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5 It is to be understood that the above-described embodiments are illustrative of only a 

few of the many possible specific embodiments which can represent applications of the 
principles of the invention. Numerous and varied other arrangements can be readily devised 
accordance with these principles by those skilled in the art without departing from the spirit 
and scope of the invention. 



in 
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We claim : 

1 A method for mobile station location estimation comprising the steps of: 

a. obtaining range measurements between said mobile station and a base station; 

b. identifying whether said base station is line of sight with said mobile station or 
non-line of sight with said mobile station at the time at which a mobile location estimate is made; 

c. correcting non-line of sight range measurements for a base station identified as 
non-line of sight with said mobile station in step b to determine reconstructed line of sight range 
measurements; 

d. repeating steps a through c for a predetermined number of said base stations; 

and 

e. determining said mobile station location estimation from said reconstructed line 
of sight range measurements determined in step c or said range measurements determined in step 
a for an identified line of sight base station in step b, or the combination of said reconstructed line 
of sight range measurements determined in step c, and said range measurements determined in 
step a for an identified line of sight base station in step b. 



2. The method of claim 1 wherein step b comprises the steps of: 
obtaining line of sight range measurements between said mobile station and a base 
station without noise, 

obtaining noisy line of sight range measurements between said mobile station and 
said base station, 
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predetermining a first standard deviation of the difference of said line of sight 
range measurements with said noisy line of sight range measurements, 

smoothing said range measurements determined in step a; 

determining a second standard deviation of the difference between said smoothed 
range measurements and said noisy line of sight range measurements; and 

discriminating between said base station being line of sight or said base station 
being non-line of sight from said first standard deviation and said second standard deviation, 
wherein said base station is determined to be non-line of sight when said second standard 
deviation is greater than said first standard deviation and line of sight when said second standard 
deviation is on the order of said first standard deviation. 

3. The method of claim 2 wherein the range measurement obtained in step a is 

represented by: 

r m (t, ) = £„(/,) + n m (/, ) + NLOS„ (t, ) 
for m = 1, M i = 0, ... K-l, wherein 

U (/, ) is the LOS distance between a mobile station and the m* base station in two dimensions 
which is given by: 

U (/,) = \x(t, ) +j * y( tl ) - Xm -j *y m \- 
J = yf-J , | | is absolute value, 

*('. ), X'. ) and (x„, ,y m ) are respectively the coordinates of the mobile station at time, /, , and 
those of the m' h base station; n m ( tj ) represents conventional measurement noise such as additive 
white Gaussian measurement noise and NLOS m (/,) represents NLOS measurement error at time 
/,.; and M is the total number of base stations; and K is the total number of time samples. 
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4. The method of claim 3 where the range measurement is smoothed by modeling: 

r m (td = 2>. («)'," 
and solving for the unknown coefficients, {a m (»)} ^ with a least squares technique. 

5. The method of claim 4 wherein the second standard deviation is represented by 



wherein 



n-0 



6. The method of claim 1 wherein step b comprises the steps of: 
estimating coordinates of said mobile station from said range measurement 
obtained in step a over time; 

calculating a range measurement from said estimated coordinates; 
determining a residual from the difference of said range measurement obtained in 
step a and said calculated range measurement; 

counting the number of times the residual is the greatest at each base station for 
each time instant; and 
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defining said base station as non-line of sight from the base station which has the 
greatest value of the number of times the greatest residual was counted. 

7. The method of claim 6 wherein said estimated coordinates are represented by 
^CXj^COat each instance of time /, , said estimated coordinates are determined as least 



M 



squares estimates to ^ = X W. ) " A»C, ))"' 

where £.(/,) = | x{t i )-x m + J* j>(,, ) - j*y m \ . 

8. The method of claim 1 wherein step c comprises the steps of: 

determining a value of maximum noise deviation and standard deviation from said 

range measurements obtained in step a and a predetermined line of sight range measurement with 

negligible noise; 

smoothing said range measurements obtained from step a; 
graphing a curve of said smoothed range measurements; 

determining a point of maximum deviation of said range measurement below said 

curve;-. 

displacing said curve downwards to pass through said point of maximum 

deviation; and 

displacing said curve upwards by said value of said maximum noise deviation, 
thereby providing said reconstructed range measurement 
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9. The method of claim 1 wherein said steps a through c are repeated for at least 
two base stations and further comprising the step of determining angle arrival information, 

wherein said mobile station location is estimated from range measurements or 
reconstructed line of sight range measurements of said two base stations and said angle arrival 
information. 



10. The method of claim 1 wherein steps a through c are repeated for three base 

stations. 

1 1 A system for mobile station location estimation comprising: 
means for obtaining range measurements between said mobile station and a 
plurality of base stations; 

identifying means for identifying whether each of said base stations is line of sight 
with said mobile station as a line of sight base station or non-line of sight with said mobile station 
as a non-line of sight base station; 

correcting means for correcting said range measurement for each of said non-line 
of sight base stations to determine a reconstructed line of sight range measurement; and 

estimating means for determining said mobile station location estimation from said 
reconstructed line of sight range measurements or said range measurements for said line of sight 
base station, or the combination of said reconstructed line of sight range measurements, and said 
range measurements for said line of sight base stations 
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12. The system of claim 1 1 wherein said identifying means comprises: 
means for obtaining a line of sight range measurement without noise between said 
mobile station and each of said base station; 

means for obtaining a noisy line of sight range measurement between said mobile 
station and each of said base station; 

means for predetermining a first standard deviation of the difference of said line of 
sight range measurement with said noisy line of sight range measurements, 
,. means for smoothing said range measurements; 

means for determining a second standard deviation of the difference between said 
smoothed range measurements and said noisy line of sight range measurement; and 

means for discriminating each of said base stations as being line of sight or being 
non-line of sight from said first standard deviation and said second standard deviation, wherein 
said base station is determined to be non-line of sight when said second standard deviation is 
significantly greater than said first standard deviation and line of sight when said second standard 
deviation is on the order of said first standard deviation. 

13. The system of claim 12 wherein the range measurement is represented by: 
r m (/, ) = L m (/, ) + Mm (/, ) + NLOS m (j t ) 

forw= 1, ...,M / = o, ... K-l, wherein 

U (t, ) is the LOS distance between a mobile station and the m* base station in two dimensions 
which is given by: 

L m (/, ) = !*(/,. ) +j * y( ti ) -x m -j *y„\- 
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j - V-7 , I I is absolute value, 

x(ti ), y(Ji ) and (x m , are respectively the coordinates of the mobile station at time, /, , and 
those of the m th base station; n m (t; ) represents conventional measurement noise such as additive 
white Gaussian measurement noise and NLOS m (/,) represents NLOS measurement error at time 
and M is the total number of base stations, and K is the total number of time sampU 



les. 



14. The system of claim 13 where the range measurement is smoothed by 

modeling: 

U(t,) = £a m (//)/; 
and solving for the unknown coefficients, {a m (/;)} ^ with a least squares technique. 

15. The system of claim 14 wherein the second standard deviation is represented 

by 



& i=0 



wherein 



,V-/ 



n-0 



16. The system of claim 12 wherein said identifying means comprises: 
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means for estimating coordinates of said mobile station from said range 
measurements from a plurality of base stations received over time; 

means for calculating a calculated range measurement from said estimated 

coordinates; 

means for determining a residual from the difference of said range measurements 
and said calculated range measurement; 

means for counting the number of times the residual is the greatest at each base 
station for each time instant; and 

means for defining said base station as non-line of sight from the base station 
which has the greatest value of the number of times the greatest residual was counted. 

17. The system of claim 14 wherein said estimated coordinates are represented by 
*«(',)..?«(',) at each instance of time t, , said estimated coordinates are determined as least 

squares estimates to F t =£(r.(/,) - L m {t,)) 2 

where L m (t, ) = | f (/j ) - Xm + j* y (i )-j* y J. 

18. The system of claim 12 wherein said estimating means comprises: 

means for determining a value of maximum noise deviation and standard deviation 
for each of said range measurements and a predetermined line of sight range measurement with 
negligible noise; 

means for smoothing said range measurements; 
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means for graphing a curve of said smoothed range measurements; 
means for determining a point of maximum deviation of said range measurements 
below said curve; 

means for displacing said curve downwards to pass through said point of 
maximum deviation; and 

means for displacing said curve upwards by said value of said maximum noise 
deviation, thereby providing said reconstructed range measurement. 

19. The system of claim 12 further comprises means for obtaining angle arrival 
information wherein said mobile station location is estimated from range measurements or 
reconstructed line of sight range measurement of said base stations and said angle arrival 
information. 
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